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ABSTRACT 

Strong gravitational lensing provides some of the deepest views of the Universe, enabling 
studies of high-redshift galaxies only possible with next-generation facilities without the lens¬ 
ing phenomenon. To date, 21 cm radio emission from neutral hydrogen has only been detected 
directly out to z ~ 0.2, limited by the sensitivity and instantaneous bandwidth of current ra¬ 
dio telescopes. We discuss how current and future radio interferometers such as the Square 
Kilometre Array (SKA) will detect lensed HI emission in individual galaxies at high redshift. 
Our calculations rely on a semi-analytic galaxy simulation with realistic HI disks (by size, 
density profile and rotation), in a cosmological context, combined with general relativistic ray 
tracing. Wide-held, blind HI surveys with the SKA are predicted to be efficient at discovering 
lensed HI systems, increasingly so at z > 2. This will be enabled by the combination of the 
magnihcation boosts, the steepness of the HI luminosity function at the high-mass end, and 
the fact that the HI spectral line is relatively isolated in frequency. These surveys will simul¬ 
taneously provide a new technique for foreground lens selection and yield the highest redshift 
HI emission detections. More near term (and existing) cm-wave facilities will push the high 
redshift HI envelope through targeted surveys of known lenses. 
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1 INTRODUCTION 

The evolution of neutral cold gas (T < 10'* K) is fundamen¬ 
tal to our understanding of galaxy evolution. Our current cosmo¬ 
logical picture is that following recombination at z ~ 1100, the 
Universe was predominately neutral until reionisation occurred be¬ 
tween z ~ 6 — 10, driven (in poorly constrained proportion) by the 
onset of star formation and supermassive black hole accretion in the 
earliest galactic haloes. The intergalactic medium is rapidly ionized 
about 1 Gyr after the Big Bang, with the vast majority of remaining 
neutral hydrogen found in galaxies where it is sufficiently shielded 
and/or replenished (e.g. Barkana & Loeb, 2001; Fan et al. 2006). 

Within galactic haloes, cold gas is a significant mass compo¬ 
nent and evolutionary ingredient: it is the pristine fuel for star for¬ 
mation, as well as the exhaust of supernova explosions. The biggest 
obstacle in uncovering the precise evolutionary role of HI is that 
direct HI observations are currently restricted to the modern Uni¬ 
verse (z < 0.2). Many important questions remain observationally 
unanswered as a result: does the cold gas continually condense onto 
galactic disks, resupplying that which is driven out by AGN and su- 
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pernovae (e.g. Dave et al. 2013, Morganti et al. 2013)? How do HI 
morphology and dispersion evolve out to cosmic noon (z ~ 2) and 
beyond? What is the link between the evolution of HI and cosmic 
star formation history of the Universe (e.g. Obreschkow & Rawl¬ 
ings 2009)? Addressing these questions demands HI observations 
at look-back times much greater than ~ 10 dynamical times (i.e., 
at z > 0.2). Since Hi extends to large galactocentric radii, high-z 
HI observations may also enable studies of the halo structure evo¬ 
lution and allow an important measurement of the baryon angular 
momentum (Obreschkow et al. 2015). 

Over the next decade we are set to revolutionize our view of 
the cosmic evolution of neutral hydrogen. A broad range of radio 
telescopes are being developed to tune into the rest-frame 21 cm 
spectral line out to large redshifts. There are a number of upcom¬ 
ing radio facilities, including the Australian Square Kilometre Ar¬ 
ray Pathfinder (ASKAP; Johnston et al., 2007), MeerKAT (Jonas, 
2009), SKA (e.g. Dewdney et al. 2013), as well the upgraded Karl 
G. Jansky Very Large Array (VLA), Westerbork Synthesis Ra¬ 
dio Telescope (WSRT) with the APERture Tile In Focus project 
(APERTIF, Verheijen et al. 2008), and the Giant Metrewave Radio 
Telescope (GMRT, Ananthakrishnan, 2005). These telescopes will 
all perform deep HI surveys, some of which will detect individual 
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galaxies out to redshifts of 2 ~ 1 and beyond. These telescopes 
will enable the study of HI in galaxies with sample sizes and at 
distances not practically possible with current facilities. This will 
renew our understanding of neutral hydrogen and its role in galaxy 
evolution. While HI absorption studies have made important con¬ 
tributions to this field (e.g. Gupta 2006, Allison et al. 2015), there 
remain significant sources of uncertainty in absorption systems that 
limit extrapolation to the cosmic evolution of HI in galaxies (e.g. 
spin temperature, gas morphology). 

One of the byproducts of the significantly larger survey speed 
and extended redshift range of next-generation instruments will be 
the dramatically increased probability that an H I-detected source is 
strongly lensed by a foreground galaxy. Here we use the canonical 
definition of strongly lensed as a total magnification p > 2. This 
corresponds to a factor jd > i saving in integration time to detect a 
given source, provided it is unresolved in the image plane. This will 
enable the detection of HI at even greater cosmological distances 
than otherwise expected with future instruments. 

In this letter, we use an N-body simulation with semi-analytic 
prescriptions for the spatially-extended HI sources to predict the 
magnification statistics and resultant number counts of strongly 
lensed Hi emission. Our simulation uses a 150 deg^ field of view 
mock observing cone out to a redshift limit of Zaource < 4, en¬ 
abling the prediction of detection rates of future HI surveys. We as¬ 
sume cosmological parameter values of Hm = 0.25, Ha = 0.75, and 
Ho = 73 kms“^Mpc~^ (Millennium run cosmology, Springel et 
al. 2005). The best estimate of the cosmological parameters today 
(e.g. the latest Planck cosmological parameters, Planck collabora¬ 
tion 2015) are slightly different, but this has a marginal difference 
on both the semi-analytic results (Guo et al., 2013) and magnifica¬ 
tion estimates, and thus on our predictions. 


2 SIMULATION 
2.1 Mock sky 

The traditional approach in estimating the lensing probability (also 
known as the lensing opacity) is to use analytic volume densities 
of both the foreground and background source populations and in¬ 
tegrate for a given redshift range. This approach has been readily 
applied in lensing predictions for point-like sources (e.g. optical 
quasars, high redshift supemovae, Oguri & Marshall 2010; com¬ 
pact radio AGN, McKean et al. 2015), however, the large spatial 
extent and non-trivial density-velocity structure of HI disks pose 
a challenge to estimating total magnification for HI sources. This 
challenge is further complicated by the requirement of modelling 
large cosmological volumes (^ 10® galaxies) in order to obtain 
sufficiently converged statistics to support actual galaxy surveys. 
We address this challenge using the S®-SAX model (Obreschkow 
et al. 2009), a semi-analytic model (SAM) extending the SAMs 
by Croton et al. (2006) and De Lucia & Blaizot (2007), relying on 
the Millennium W-body simulation of cosmic structure (Springel et 
al., 2005). S®-SAX is the only model to date that provides a mock 
sky of millions of HI disks with realistic spatial distributions and 
velocity profiles (Obreschkow et al., 2009b). This model is well- 
suited to generate HI lensing predictions since the mock galaxies 
have all the properties needed for both the source and lens galaxies. 
Explicitly, (i) the mock galaxies are correctly clustered due to the 
underlying Millennium simulation; (ii) the HI (needed in source- 
galaxies) has realistic spatial distributions, accounting for the large 
HI disk size (> stellar effective radius) and central troughs; (iii) the 
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Figure 1. Top: Image and source plane channel maps of a 2 : = 1.86, 
Mhi = 2.7 X 10® Mq H I disk (frequency resolution = 200 kHz), demon¬ 
strating how the approaching HI emission overlaps with the inner lens caus¬ 
tic and results in a full HI Einstein ring in some channels. The lens is at a 
redshift of 2] = 0.58, with a total mass Muaio = 3.6 X lO'^® Mq and bulge 
ellipticity en = 0.24. Each frame is 6.4x6.4 arcsec®. The caustics/critical 
curves are plotted in red. Bottom: Lensed and unlensed HI profiles of the 
same system at higher frequency resolution (50 kHz). Significant differen¬ 
tial magnification is seen between channels to the extent that the probability 
of detection is higher for a narrow (~50 kHz) channel width, rather than the 
full HI line width. This argues that the detection statistics presented in this 
paper are likely to be conservative if source detection algorithms are also 
run on narrow channel widths. 


HI has realistic velocity profiles that also satisfy the Tully-Fisher 
relation and local Hi velocity function (Obreschkow et al. 2013); 
(iv) the galaxies have multiple components (halo-l-disk-l-bulge) with 
realistic mass distributions (needed for lens galaxies), (v) the prop¬ 
erties of the mock galaxies evolve with redshift according to the 
prescriptions of Croton et al. (2006) and Obreschkow et al. (2009). 
A velocity dispersion of CTv.hi = 10 kms“^ is added to the SAM¬ 
generated HI velocity profile, consistent with observed properties 
of nearby HI galaxies (e.g. de Blok et al. 2008). The typical HI ve¬ 
locity dispersions and scale heights potentially increase with red¬ 
shift, similarly to the observed increase in the molecular velocity 
dispersion (e.g. Forster Schreiber et al. 2006). We decided to ne¬ 
glect this effect, since the rest-frame channel widths and pixel sizes 
of our lensing calculations are generally a few times larger than 
the dispersion scale and thickness of the disk. Crucially, the use 
of an N-body simulation and full source-to-lens plane ray-tracing 
enables the HI flux of all lensed galaxies to be split into individ¬ 
ual frequency channels, resulting in a realistic prediction of the ob¬ 
servables (e.g. peak channel flux, velocity line widths) and hence 
optimal detection strategies. 


2.2 Lens candidate selection 

Based on the S®-SAX model we produce five statistically indepen¬ 
dent mock observing cones (as described in Obreschkow et al., 
2009b), each covering a field-of-view of 30 deg® and extending 
from 2 = 0 to 2 = 4 (30 deg® corresponds to the cross-section 
of the Millennium box at 2 ~ 4). Together these five cones rep¬ 
resent a mock sky covering 150 deg® with negligible cosmic vari¬ 
ance. From this sky we extract all lens-source pair candidates with a 
background source integrated HI flux 5^*3 > 100 )ijy km s“^ and 
foreground halo virial mass Mvir,i > 5 x 10^^ Mq. In addition, we 
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only consider lens candidates with a source redshift ^source < 4 
and an intrinsic lens-source impact factor < (20 arcsec) AND 
(0.15 Rviv), where R^ii is the foreground halo virial radius. These 
cuts dramatically decrease the computational demand and have a 
negligible effect on the inferred lensing statistics (following a trial 
and error process). For example, the majority of lensing galaxies 
have Einstein radii (and hence deflection angles) of order a few 
arcseconds (e.g. Treu 2010, and references therein). The above se¬ 
lection criteria result in ~ 10® lens-source pair candidates, i.e. ap¬ 
proximately 2 per arcmin^. 

2.3 HI magnification factor 

The lens models are split into three components: (a) a dark matter 
halo, (b) bulge, and (c) a disk component. The dark matter halo 
is assumed to have a spherically-symmetric NFW density profile 
with a concentration parameter derived from a redshift and halo 
mass dependent fit (Hennawi et al. 2007), the distribution of which 
is consistent with empirical results from X-ray and strong lensing 
observations of galaxy clusters (Comerford & Natarajan 2007). The 
latter consistency is particularly important since magnification bias 
will preferentially select high surface mass density foreground dark 
matter haloes. 

The bulge and disk effective radii are derived directly from 
the self-consistent dynamics of the semi-analytic prescription to all 
galaxies in the simulation. The bulge is set to a de Vaucoleurs den¬ 
sity profile with an ellipticity drawn from a Gaussian distribution 
with mean and Icr standard deviation of ~ 0.3 ± 0.3 (truncated 
between 0 < pe < 0.98), consistent with observational proper¬ 
ties of bulges and previous approaches (Jorgensen & Franx 1994, 
Mdller et al. 2007, Oguri & Marshall 2010). The halo-l-bulge den¬ 
sity profiles used in SAX-Lens are realistic not only in that the 
distribution of their concentration parameters match available ob¬ 
servations; but also in that global Hi linewidths (outer rotation) 
and CO linewidths (inner rotation) are both consistent with ob¬ 
servations (Obreschkow et al. 2009b). The disk is modeled as a 
spherically-symmetric component with an exponential density pro¬ 
file. It is sub-dominant at all radii and adds a marginal fraction of 
total convergence to the system as seen in Moller et al. (2007). We 
do not consider the additional convergence contributed by satel¬ 
lite or neighbouring galaxies associated with the foreground lens, 
which may add up to several tens of percent to some of the largest 
image separations (and hence magnification, Oguri 2006), however 
this effect decreases for lower mass, lower redshift lenses as dis¬ 
cussed by Moller et al. (2007). 

With both source and lens profiles defined, we perform the 
general-relativistic ray-tracing using the glafic software pack¬ 
age (Oguri 2010). To this end, the HI emission of all sources is split 
into 50 kHz channels (intrinsic velocity widths of ~ 11—53 km s“^ 
for 2 = 0 — 4), each of which is individually lensed. We thus have 
noise-free image plane (lensed) and source plane (unlensed) ob¬ 
serving cubes (two spatial dimensions, one frequency dimension) 
of all HI sources. Fig. [T] shows example source and image plane 
channel maps of a 2 ~ 2 source (at 200 kHz channel resolution 
for clarity). The associated lensed and unlensed global profiles are 
shown in the bottom panel. Optimal detection strategies for such a 
distorted profile are discussed in Sec. |3.3| The global magnification 
factor /i is defined as the ratio between the integrated HI fluxes of 
the magnified image and the unmagnified source. All systems with 
/r > 2 are written to a catalogue, referred to as the Semi-Analytic 
eXtragalalactic Lens (SAX-Lens) catalogue. This catalog contains 
45, 734 lens-source pairs, i.e. about 300 per deg^. 



Figure 2. Top: Redshift distribution of all foreground lens and background 
source pairs with total HI magnification /r > 2 (ff units: deg“^). Bottom: 
The fraction of HI sources that are predicted to be strongly lensed (lensing 
opacity), given a velocity-integrated flux threshold (see legend). The grey 
shadings indicate the Ic Poisson error, given the full 150 deg^ simulation. 
This forecasts the efficiency of future HI surveys at selecting lensed sys¬ 
tems, paificularly if cross-correlated with massive foreground galaxies at 
lower redshift. 


3 PREDICTIONS 

3.1 Survey-independent lens statistics 

The primary output of the methodology laid out in Sec. is the 
SAX-Lens catalogue of simulated HI sources that are magnified 
by a factor of 2 or more. In Fig.|^(top panel) we show the redshift 
distribution of both the foreground lenses and the background HI 
sources. This demonstrates the lensing probability only becomes 
appreciable above 2 > 0.5, well above the highest redshift Hi 
emission detection to date (z ~ 0.2, Catinella et al. 2008). The dis¬ 
tribution of strongly lensed HI disks with Shi.s > 100 ply km s“^ 
peaks at 2 ~ 1.5—1.8, which corresponds to an observed frequency 
ofiZobs ~ 570—510MHz. The bottom panel in Fig.j^demonstrates 
that higher redshift sources are, on average, more magnified due to 
higher alignment probability and smaller apparent source size. This 
means that surveys designed to detect these sources will be highly 
biased towards strongly lensed systems. It is likely that future ra¬ 
dio telescopes able to detect lensed Hi sources, will, due to the 
wide (> 2 : 1) instantaneous band widths of modern receivers, 
simultaneously detect the HI emission or absorption of the fore¬ 
ground lens. This will greatly expedite the followup observations 
of lensed sources discovered with an entirely different selection 
function to other lens catalogues (e.g. CLASS, Myers et al. 2003; 
SLAGS, Bolton et al. 2006). 

In Fig. we show some of the statistical properties of the 
lens-source pairs, including the foreground lens halo (total) mass, 
impact parameter (angular distance between lens and source cen¬ 
troids), magnification (colour scale), as well as the HI disk angular 
size. The ellipse axis ratios correspond to that of the background 
Hi source (with min=0.01). The ellipse position angle represents 
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Figure 3. Total magnification as a function of foreground lens virial mass 
and impact parameter. The magnification is clipped to /r = 10 and only 
every 10**' lensed source in the full 150 deg^ simulation is shown. Ellipses 
show the intrinsic angular size of the background HI disks, with a I ai'csec 
reference shown in black (below inset). This shows that the highest magni¬ 
fication sources tend to have small angular size. The ellipse position angles 
represent the angle between the HI semi-major axis and the impact vector. 
Vertical (horizontal) semi-major axes represent aligned (perpendicular) im¬ 
pact vectors and HI semi-major axes. The dark grey line shows the fit to the 
upper envelope (see text). Inset: Total magnification distribution of /r > 2 
sources based on full 150 deg^ simulation (bin width = 2). 

the angle between the impact vector and HI semi-major axis (see 
caption). Fig. shows a rough envelope above which virtually no 
sources are strongly lensed (marginalized over a number of param¬ 
eters, e.g. redshifts, halo concentration parameter). An approximate 
fit to this envelope is: impact ~ 0.50 + (Miens/10^^ ®® M©)** ®® 
arcsec. 

The highest magnification sources (/r > 10) in the SAX-Lens 
catalogue have an Einstein ring (see Fig.[TJ in at least one of their 
channel maps and thus have the largest apparent solid angles. How¬ 
ever, current and future radio facilities will not have sufficient angu¬ 
lar resolution (for the vast majority of their collecting area) at the 
appropriate frequency to spatially resolve the HI Einstein rings. 
This maximises their probability of detection since the surveys will 
benefit from the total Hi magnification boost. This also implies 
that macro lens models derived from optical/NIR/mm observations 
will be required to constrain the intrinsic HI disk mass and size. 
Of all the lensed HI disks in the SAX-Lens catalogue, just 6 per¬ 
cent have an extent larger than 0.5 arcsec. This implies that even 
SKAl-MID, with a < 0.5 arcsec PSF for Uabs ^ 1 GHz, will not 
resolve the vast majority of lensed HI disks. Of course, the sources 
that are resolved will be of great interest for individual followup to 
perform more detailed modelling of the astrophysical properties. In 
this paper we limit the scope to the detection statistics, assuming 
all sources are unresolved. 


3.2 Survey-specific predictions 

Having presented the basic statistical properties of the lensed HI 
systems, we now turn to predicted detection rates with planned sur¬ 
veys. There are a wide range of HI surveys planned over the course 
of the next decade which will dramatically expand our view and 
understanding of HI in galaxies. In Tab. [T] we list the predicted 
number of HI lenses (and mean redshift) that will be detectable for 
some of the major future surveys. We define the detection thresh¬ 
old as .Shi > 5a, where .Shi the mean HI profile flux in units of ly 
and the noise corresponds to that expected in a channel width equal 


Table 1. Predictions for the number of Hl-detectable lensed systems in 
future surveys. As described in the main text, a detection is defined as 
S'mean > 5cr. For the single pointing surveys (CHILES, LADUMA), the 
area (H) corresponds to the 1.4 GHz field of view, which scales as (1 -|- z)'^. 
SKA surveys assume the SKAl-MID specifications, ~1 year integrations 
and indicative survey areas. However, we note that SKAI-LOW could cover 
the z > 2 space if its upper frequency limit were to be increased to 
~470 MHz. The SKAl-MID specifications are drawn from Dewdney et 
al. (2013), but the total number of antennas decreased from 254 to 197 in 
line with the 2015 re-baselining process. The quoted SKAl-MID survey 
sensitivities are assumed for 1420.4 MHz and 1050 MHz, where the lat¬ 
ter is the highest frequency in Band 1. Sensitivities for both SKA Band 1 
and 2 scale as (1 -|- 5:)“* due to the increase field-of-view but fixed survey 
area. Multi-wavelength cross-coiTelation may reveal many more HI lenses, 
as suggested by the 3a predictions in parentheses. Note that the APERTIF 
surveys resulted in zero detections in our simulation, so detection probabil¬ 
ity estimates were not possible. For reference, the approximate number of 
strongly lensed systems known at present is ~ 500. 

“ COSMOS HI Large Extragalactic Survey (Fernandez et al., 2015) 

^ Looking at the Distant Universe with the MeerKAT Array (Holwerda. Blyth & Baker, 2012) 

^ Deep Investigations of Neutral Gas Origins (Meyer, 2009) 


to the HI profile FWHM for each individual galaxy. The latter are 
determined from the 100 kHz channel sensitivities listed in Tab.[T] 
A major caveat in any such prediction is the ability to distinguish 
lensed from non-lensed sources, given that they will be unresolved 
(in HI ) and the vast majority will only be lensed by a factor of 
two (see inset, Fig.|^. Multi-wavelength data will therefore be es¬ 
sential to unambiguously identify low (fi < 4) magnification HI 
lenses, but may also be used to lower the detection threshold and 
discover many more lensed systems, as suggested by the number of 
3a detections indicated in parentheses. 

Table [T] shows that none of the pre-SKA surveys are particu¬ 
larly optimized for detecting large numbers of lensed HI systems 
(e.g. ~ 4 — 20 detections for DINGO and LADUMA). However, 
lensed HI systems will increase the number counts and enhance 
these already planned HI surveys beyond their original science 
cases. The predicted mean redshifts of these lensed HI sources 
are much higher than the current record of z ~ 0.2 (except for 
DINGO-D). This combined with the magnification boosts implies 
that it will be possible to detect high redshift systems with much 
lower HI mass than that possible for non-lensed systems. Further¬ 
more, the relative isolation of the HI line, the steep high-end HI 
mass function and the systematically larger magnification factors 
at higher redshift suggest that some of the most distant HI sources 
will be detected and identified within a fraction of the total survey 
duration. See Serjeant (2014) for a discussion on using the high- 
mass end of the HI mass function to select lenses. The above fac¬ 
tors, combined with the increasing HI lens opacity at higher red¬ 
shift (see Fig. bottom panel) will make a significant impact on 
Hhi measurements at 2 > 2. Indeed, even 2 ~ 1 Hhi values with 
pre-SKA facilities may require strong-lensing correction factors. 
Tab.[T]also suggests that targeted surveys for lensed HI may be the 
optimal strategy in the pre-SKA era. 
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The SKAl surveys will revolutionize our view of HI in galax¬ 
ies at high redshift. The predictions suggest ~ 0.8 — 2.2 x 10"* 
lensed detections, well over an order of magnitude larger than the 
presently known lens population at any wavelength. The mean and 
maximum redshifts of sources detected in these surveys will pro¬ 
vide important insights on the evolution of cold gas towards cosmic 
noon (2 ~ 2), a perspective not possible with any other facility. 

The full SKA is expected to detect up to ~ 10® non-lensed 
HI galaxies out to 2 ~ 1 through surveys aimed primarily at high 
precision measurements of baryonic acoustic oscillations (Rawl¬ 
ings et al. 2004, Maartens et al. 2015, Bull et al. 2015). However, 
even with the SKA, strong gravitational lensing will be required to 
characterize Milky Way like galaxies to the epoch of peak cosmic 
star formation and AGN activity at 2 ~ 2 — 3. Stacking methods 
have been proposed to pursue this goal, however, the uncertainties 
and associated contamination using these methods are unclear. Fur¬ 
thermore, calibration and source modelling limitations may place 
upper limits on the sensitivity of stacking methods (Makhathini 
et al., in prep.). Well-modeled strong gravitational lenses, with 
the benefit of excellent multi-wavelength coverage from ALMA, 
LSST, JWST and the E-ELT will clearly provide an excellent op¬ 
portunity to push back this observational frontier. While these sur¬ 
veys will dramatically enhance our knowledge of high redshift cold 
gas and disk formation, a further valuable contribution will be a 
new probe of the foreground lens population using a unique se¬ 
lection function that is not biased by dust attenuation (in both the 
source and lens galaxies). 


3.3 Optimal detection strategy 

In all of the lensed statistics presented above, the detection thresh¬ 
olds assume a channel width matched to the FWHM of the HI line 
width. This optimizes the S/N ratio and hence maximizes detection 
probability for the majority of lensed systems. However, there are 
cases when a narrow channel width has a significantly larger mag¬ 
nification than that of the total intensity map. This is simply due to 
the channel flux solid angle and position relative to a caustic. Such 
differential magnification is seen in the relative magnification of 
the CO (1-0) channel maps of a 2 ~ 2.3 lensed disk galaxy (Deane 
et al. 2013a), as well as the order of magnitude difference in mag¬ 
nification of the radio core to CO emission (Deane et al. 2013b). 
Differential magnification can be exploited to increase the detec¬ 
tion probability by running spectral line source finders with narrow 
channel widths as well as those optimized to the full galaxy line 
width. An illustration of this effect is seen in Fig.[T]which shows a 
2 ~ 2 Hi source with total magnification /i ~ 3.3, however one 
particular 50 kHz channel is magnified by a factor n > 10. An 
appropriate channel selection yields a higher probability of detec¬ 
tion than matching to the line width. Clearly, there is a great deal 
of optimization to be explored in both targeted and blind lensed HI 
surveys, which we investigate in an upcoming paper. 


4 CONCLUSIONS 

The use of an N-body simulation with applied semi-analytics en¬ 
ables the following predictions relevant to lensed Hi. 

(i) A full ray-tracing analysis of 10® 2 < 4 lens candidates in the 

5 X 30 deg® simulation reveal ~46,000 strong lenses (/r > 2) with 
unlensed velocity-integrated Hi fluxes S'hi > 100 pjykms“^. 


(ii) The predicted number counts imply that the dramatic in¬ 
crease in sensitivity and frequency coverage of upcoming radio fa¬ 
cilities will enable HI detection at redshifts where strong lensing 
has been highly beneficial to galaxy evolution studies that utilise 
other spectral lines (e.g. CO, [Cll], [Olll]) as well as continuum 
observations. 

(iii) Planned HI surveys (LADUMA, DINGO) on future tele¬ 
scopes will detect some (~ 4 — 20) lenses, which in the case of 
LADUMA may require consideration when determining Ohi( 2 : ~ 
1). However, none of these surveys are particularly optimal for de¬ 
tecting lensed systems. Lensing does however enable the detection 
much lower HI mass systems at high redshift, enhancing the scien¬ 
tific output of planned surveys. 

(iv) Current telescopes (e.g. GMRT, VLA) fitted with their new 
broadband receivers, and most definitely future telescopes (e.g. 
MeerKAT), are very well placed to detect lensed Hi in targeted 
surveys and expand our view of cold neutral gas to higher redshift. 

(v) The depth and area of the SKAl-MID surveys will lead to 
efficient lens selection at high redshift. Overall, these surveys are 
predicted to detect over an order of magnitude more Hi lenses 
(~ 0.8 — 2.2 X 10'*) than the presently known lens population at 
any wavelength. 

(vi) The large solid angles covered by Hi disks leads to sig¬ 
nificant differential magnification in some lensed sources. Appro¬ 
priately optimized three-dimensional source finders may therefore 
enhance the detection rates presented here. 
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